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ABSTRACT 


FANTASTIC is an acronym for Failure Analysis Nonlinear Thermal 
and Structural Integrated Code. This nroqram has beer) anvolonp') 
by Failure Analysis Associates, Palo Alto, California for MSFC in 
order to improve the accuracy of solid rocket motor nozle 
analysis. Its first version has boon released and received by 
the Thermal Systems Branch of MSFC. FANTASTIC has three 
modules : 


a . 

FACT: 

It 

is 

the 

thermochern ica 1 analysis modul 

h. 

FAHT : 

It 

is 

the 

heat transfer analysis module 

c . 

FAST: 

It 

is 

the 

structural analysis module 


All modules have keywords for input of data. Work is in proqress 
t 01 the ver if ic at ion of F AHT modu 1 o • Th is is b<* i nq don© by us j. nq 
data for various Problems with known solutions as inputs to the 
bAHT modulo. The information obtained by runninq the.se proqrams 
is used to ninnoint Problem areas of this code and passed on to 
the developer for removinq huqs from this code. As a result of 
this Procedure , failure Analysis Associates have revised the 
first version of the FANTASTIC code and a now improved version 
has been released and received by the Thermal Systems Branch. 

Fur. thor work to verify the new version is continued: 


i 

i 

i 


I 
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1 . Introduct ion 

Failure Analsysis Associates, Palo Alto, California has developed 
an advanced computer code for USFC to improve the accuracy of the 
solid rocket motor analysis. The commuter code is called 
"FANTASTIC", an acronym for Failure Analysis Nonlinear Thermal 
and Structural Integrated Code. The code incorporates a higher 
order inteqration scheme for model inq steep temperature and 
strain Gradients. It has a modular design and is intended as a 
multi-purpose thermostructural analysis code. 

FANTASTIC has not been verified before. In order to us a this 
code for the intended purpose outlined above, it is necessary to 
verify it. The purpose of this report is to give an overview of 
the procedure followed to verify this code. 

2 . Procedure 

FANTASTIC consists of three modules: 


1 . 

FACT: 

It 

is 

the 

thermochem ica 1 analysis module 

2. 

FAHT: 

It 

is 

the 

heat transfer analysis module 

3. 

FAST : 

It 

is 

the 

structural anal vs is modul e 


Because FANTASTIC can be used for a stand-alone thermal or 
structural analysis, FAMT module was chosen for verification 
purpose. FAHT module accents data under numerous keywords. All 
keywords of FAHT are divided into three blocks: 

1. PARA'-'ETFR Block: This block is used to nrovide a title 

and analysis procedure (i..e., steady or transient etc.) options. 

2. MODEL Block: This block is used to provide the details 

of the model geometry, material, and boundary conditions. 

3. INCREMENT Block: This block is used to provide the 
loading history in the form of various types of incremental 
1 oads . 

A list of the keywords under each block is given in Table 1. The 
data must be input as required by the orogram format. 
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2-1 MODEL Geometry 

The first requirement in setting un the numerical ^robl^n for 
execution by FAHT is to define the domain of solution. Thnn this 
domain is divided into apnropr iato number of "finite elements". 
There are various types of elements and an element can have nodes 
ranging from 2 (link element) to 20 ( isoparametric hexahdron 

element). A list of the names for various tynes of elements 
covered by FAHT is qiven in Table 2. As a qeneral. rule more 
accuracy is obtained by select inq eiqht noded olanar element than 
four noded nlanar element for the same curved geometry. In three 
dimensional nroblems un to 20 noded elements mav be chosen 
depend inq unon the oeometry. Each node is fixed by a set of 
coordinates which are input under an annrooriate keyword in the 
MODEL block . 

2.2 Material Description 

There are a total of 12 FAHT material, types. Each material type 
is defined by a set of properties. For example Tyne 1 - Thermal: 
Linear Isotropic material requires three properties namely mass 
density, specific heat, and thermal conductivity whereas Type 41 
— Isothropic three— component charting and decomoos inq material 
requ ires 32 properties to be specified. A list of FAHT material 
types is qiven in Table 3. 

2.3 Loading Conditions 

Loading conditions treated by FAHT include convection, enclosure 
radiation, lumped heat capacity, and heat sink/source. The 
convection option is implemented through the keyword 
*HEATCONVECTIONBC group. The radiant heat exchange option is 
implemented through the * f) NC LOS U P E RA n I AT I ON and *VIEWFACTOR 
groups. Lumped heat capacity option is implemented via the 
*LUMPEDHEATCAPACITY group. A volumetric heat sink or source can 
be defined as a function of time or temperature via the 
♦HEATGENERATION GROUP. 

2.4 Boundary Conditions 

Boundary conditions f or the boat transfer analysis may be 
specified in terms of prescribed tempers ture or heat flux along 
external boundaries or internal surfaces (e.q., contact 
r-sostance or heat flux across gaps between Riir f ac»s). This 
option 'is implemented via the *TEMPER ATURF.BC and *HF.ATFLUXPC 

groups . 
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2.5 Initial Conditions 

Initial temperatures ma y he snr>c if i ed in "AHT at each node. 
Internal heat qeneration rates can also be specif ied as functions 
of time or temnerature. 

3. Illustrative Example 


Consider a two-dimensional steady state heat transfer in a 
material with constant thermal conductivity. A schematic showinn 
the boundary conditions is q iven in Fiqure 1* ’The qeometrv has 
32 elements. Each element has four nodes. Thus the total number 
of nodes is 45. The element type from Table 2 is type number 
21, planar quadrilateral. The material type is assumed to be 
type number 1, thermal: linear isotroic. Therefore, three 

properties namely mass density, snecific heat, and thermal, 
conductivity are required to be specified. This information is 
qiven in the material library section of the FAHT user's manual. 
The input for the imol emen tat ion of FHAT for this Problem under 
suitable keywords of the PARAMETER block, MODEL block r and the 
INCREMENT block respectively is niven below. Explanations for 
.the input data are also qiven: 

PARAMETER BLOCK 

♦TITLE 

Two Dimensional Steady State float Transfer 
♦TEMPERATURE 

This keyword is used because the h^at transfer analysis is to do 
performed, in the next time period. 

*i;ndparametep. 

F'r is keyword is use"* to define the end of the PARAMFTFR block. 
MODEL Block 


♦MODE 

This qrouo is used to define nodal coordinates. A node setname 
must also be defined ^or a qrouo of nodes. 
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More NODESET is the name of the node set consistinq n f all <15 
nodes . 

1 0 . 0 0 . 0 

The first number 1, denotes a node which is assiqned number 1 
with coordinates x = 0.0, y = 0.0. The next two numbers 0.0 and 
0.0 denote the x and v coordinates of node 1. 

9 2.0 0.0 1 

The first number 9, denotes a node wh ich is assiqned number 9 
with coordinates x = 2.0, y = 0.0. The next two numbers 2.0 and 
0.0 denote the x and y coordinates of node 9. The last number 1, 
denotes the increment for nodes bwetween nodes 1 and 9. The 
coordinates of nodes 2 through 9 are oenerated automatically by 
the build in generator in the code. Ue can similarly input the 
coordinates of th^ remaining nodes as given ho low. 


■Id 

0.0 

0.25 


18 

2.0 

0.2 5 

l 

19 

0.0 

0.50 


27 

2.0 

0. 50 

1 

28 

0.0 

0.75 


36 

2.0 

0.75 

1 

37 

0.0 

1.00 


45 

2.0 

1.0 0 

1 


* ELEMENT 

This qroup is used to define element connectivity. Each element 
must be assiqned a unique element number. An element set. name 
must also be defined for each qroup of elements . 

ELSKT 21 1 

The name of the element set consistinq of all 32 elements is 
ELSET. The number 21 denotes the element tyoe number and the 
next number 1 denotes the mat^iral number for the elements in 
this set. 

1 12 11 19 7 1 


XIII-4 


ORIGINAL page is 

OK POOR QUALITY: 


The first number 1 denotes an element o f the set. The next four 
numbers 1, 2, 11, and 10 denote that the element is connected hy 

these four nodes. The next two numbers 7 and 1 denote fhat there 
are 7 similar elements 'with increment 1. This data line, 
therefore, covers inut for eight elements 1 through 3. The data 


for 

rema in i ng 

elements 

is 

i nput 

similarly as below: 

9 

10 

11 

20 

19 

7 

1 

17 

19 

20 

29 

28 

7 

1 

25 

28 

29 

38 

37 

7 

1 


*SET 


This group is used to define additional node and element sets. 
Each set must be assigned a unigue set name. A node set or an 
element set may he defined as a list of (node or element) numbers 
or as a list of (node or element) number ranges. 

NODESET 112 

here NODESET 1 is the name of a node set consisting of nodes 1 
through 9. Number 1 denotes that the set-type is node set. (If 
this number is 2, it would have denoted that the set type is 
element set). The next number 2 denotes the list tve . Here it 
means that we are going to input a range of set nodes. 

1 9 

The number 1 denotes the first node number and number 9 denotes 
the last node number (range from node 1 through 9). 

NODES ET 211 

Here NODES ET 2 is the name of the node set consisting of 
1ft, 27, 35, and 45. Number 1 denotes the sottvpo (here node a '?t) 
and the next number* 1 d o n o t e s that the list type is for a list of 
node(s) given below. 

18 27 36 4 5 

*TEMPERATUREBC 


10 

383.0 

19 

707.0 

2 8 

924 .0 

37 

10 00 . 


XIII -5 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Above we have input tor spec if ied temperatures at nodes 1 n , 1 Q > 

2 8, and 37. Another method to irout socified temperature 
boundary condition for a node c ot is shown below: 

NODESET 1 0.0 

NODESET 2 0.0 

All nodes in NODESET 1 and NODESET 2 have a specified temperature 
of 0.0. 

*PROPERT¥ 

This group is used to define material properties. 

1 1 

The first number 1 denotes the material number and the second 
number 1 denotes the material type (in this examnle - thermal: 
linear isotropic). So, three properties namely mass , dens ity , 
specific heat, and thermal conductivity are input in that order 
as below: 

1 1 1 

* END MODEL 

This keyword is used to def ine the end o f the model block. 

* PRINT 

This group is used to select the desired printed output nrintcode 
10 is used for printing of temperatures. 


NODESET 10 

With this inout temperatures for all no' q °s in the set MODESET are 
pr inted . 


*ENDINCREMENT 

This keyword is used to define the end o f the INCREMENT block 
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4 . CONCLUSIONS 

Ten problems with known solutions were usp 4 for input of data 
to FAHT. Based on the runs made and their results, it was found 
that FAHT needed revision because of buns in the program. For 
example, the program was not able to recoon ize a few kevwords for 
input of data. 

The solutions of many problems did not match the known 
solutions. This indicated that the program needed an overhaul. 
This information was massed on to the developer. Failure Analysis 
Associates. Based on these results, the developer has been abl<-> 
to remove bugs and has supplied MSC with a new version of 
FANTASTIC with workable FAHT module. 
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TABLE 1 VARIOUS DATA BLOCKS AND KEYWORDS FOR FAHT 


PARAMETER Block 

♦TITLE 

♦FACT 

'TEMPERATURE 

♦POREPRESSURE 

♦FAST 

♦TRANSIENT 

♦NONLINEAR 

♦ENDPARAMETER 

MODEL Block 

♦NODE . . . 

♦ELEMENT 

♦SET 

♦CROSSSECTION 

♦PROPERTY . 

♦MATERIALDIRECTION . 
♦TEMPERATUREBC .... 

♦PRESSUREBC 

♦HEATFLUXBC 

♦MASSFLUXBC 

♦INITIALTEMPERATURE . 

' ‘INITIALPRESSURE .... 
♦HEATCONVECTIONBC . 
♦THERMALGAPCONTACT 
♦LUMPEDHEATCAPACITY 
♦HEATGENERATION . . . 

* NONLINE ARHEATBC . . . 
♦TIMETEMPFUNCTION . . 
•ENCLOSURERADIATION . 

•VIEWFACTOR • • 

•ENDMODEL 

INCREMENT Data Block 

♦TRANSIENT 

♦NONLINEAR 

♦RESTART 

♦PRINT . . 

♦REPORT 

♦POST 

♦ENDINCREMENT 
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Table 2 : FAHT Element Types 


Type 

Element 

Number 

Number 

Description 

of Nodts 

11 

Link 

2 

21 

Planar quadrilateral 

4 

22 

Planar quadrilateral 

• 

25 

Axisymmetric quadrilateral 

4 

26 

Axisymmetric quadrilateral 

8 

31 

Solid hexahedron 

8 

32 

Solid hexahedron 

20 

41 

2-D planar interface 

4 

42 

2-D planar interface 

6 

45 

Axisymmetric interface 

4 

46 

Axisymmetric interface 

6 

51 

3-D interface 

8 

52 

3-D interface 

16 


Table 3: FART Material Types 


Type 

Number 

Material 

Description 

Number of 
Properties 

1 

2 

Thermal: linear isotropic 

Thermal: time or temperature dependent 

isotropic 

3 

6 

11 

Diffusion: linear isotropic 

4 

12 

Diffusion: time or temperature dependent 

isotropic 

8 

21 

Thermal: linear orthotropic 

8 

22 

Thermal: time or temperature dependent or- 
thotropic 

16 

31 

Diffusion: linear orthotropic 

0 

32 

Diffusion: time or temperature dependent or- 
thotropic 

18 

41 

Isotropic three-component charring and decom- 
posing materia] 

32 

42 

Isotropic time or temperature dependent three- 
component charring and decomposing material 


43 

Orthotropic three-component charring and de- 
composing materia] 


44 

Orthotropic time or temperature dependent 
three-component charring and decomposing ma- 
terial 
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♦TITLE. 


TWO-DIMENSIONAL STEADY STATE HEAT TRANSFER 
IN CONSTANT CONDUCTIVITY MEDIUM WITH INSULATED. 


AND ISOTHERMAL BOUNDARY CONDITIONS. 


IOOO 

.It 

.Jt 

INSULATED • 



-* 

* 

-* 

* 


5 37 

5 38 

5 39 

540 541 > 

542 - 

543 

5 44. ' 

545 

a 

a 


t 

t 

a 

a • 

1 

a 

1 

• 

! * 

0. 25 

924 

* 25 1 26 

* * 

T 27 
-* 

1 28 i 29 

: 30 


5 31 

-* 

5 32 

-* 

• • 

Jf. 

5 2a 

5 29 

5 30 

531 532' 

5 33 

5 34 

5 35 

1 36 



a 

a 

1 
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f 

1 

• 

1 

« 

a 
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-* 
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5 23 

5 24 
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a 

-* 

ZERO 


1 19 
I 


120 
I * 


121 
I 


• 22 


5 23 
5 
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FIGURE 1 MODEL GEOMETRY SHOWINGS NODES (*) , NODE NUMBERS 1 THROUGH 45, 
ELEMENTS 1 THROUGH 32, AND BOUNDARY CONDITIONS 
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